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optimizing tftfe at least one ligand position while allowing translation, orientation 
and rotatable bonds W the ligand to vary, and while holding the protein fixed; 
wherein said method does not involve assembly of two or more fragments to form the ligand. 
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(Amended) 4. Ttae method of claim 1, wherein said performing the pre-docking 
conformational search comprises: 

randomly generating a plurality of conformations of the ligand; 
Oi^ \, minimizing la strain of each conformation of the plurality of conformations; 

I using the strain and a solvent accessible surface area of each conformation to 

rank the conformations; and 

clustering thelconformations and retaining a desired number of top clusters of 
conformations, said retained number of top clusters of conformations comprising said 
multiple solution conformations of the ligand. 



(Amended) 6. The method of claim 5, wherein said generating the binding site image 
further comprises: 

placing a grid around the binding site of the protein; 

determining a hot spot search volume using said grid; 

determining hot spots using a grid-like search of the hot spot search volume; 

and 



for each type of hot spot, clustering the hot spots and retaining a desired 
number of top clusters of hot spots, said desired number of top clusters comprising said 
multiple hot spots to be employed by said matching. 



(Amended) 8. "Thei method of claim 7, wherein said determining the unique rigid body 
transformation comprises determining the unique rigid body transformation that minimizes: 
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position vector of a j th hot spot of the protein; 
a positfofhA^ctor of a j th atom of the at least one solution conformation; 
a 3x3 rotation mafrixs^nd 
a translation vector. 



(Amended) B. The method of claim 1 , wherein said optimizing comprises optimizing 
multiple protein\igand complex formations, said optimizing comprising: 

elmninating each ligand position having a predetermined percentage of ligand 
atoms with\a steric clash; 

ranking remaining ligand positions using an atom pairwise score with a desired 
atom score qutoff, said atom pairwise score comprising a hydrogen bonding potential 
score or a stone potential score; 

after ranking, clustering the ligand positions and selecting a top number n of 
ligand positions; and 

optimizing each ligand position of the n positions, allowing the translation, 
rotation and rotatable bonds of the ligand to vary. 

(Amended) 10. Theynethod of claim 9, wherein said optimizing comprises optimizing 
each ligand position of theVi positions using a Broyden-Fletcher-Goldfarb-Shanno (BFGS) 
optimization algorithm with said simple atom pairwise score, allowing the translation, rotation 
and rotatable bonds of the ligand to vary. 

(Amended) 11. A systemVor docking a ligand to a protein comprising: 

means for performing a pre-docking conformational search to generate multiple 
solution conformations of the ligand; 

means for generating a binding site image of the protein, said binding site image 
comprising multiple hot spots 

means for matching hot spots of the binding site image to atoms in at least one 
solution conformation of the multiple solution conformations of the ligand to obtain at 
least one ligand position relatives the protein; and 

means for optimizing the atUeast one ligand position while allowing translation, 
orientation and rotatable bonds of tfce ligand to vary, and while holding the protein 
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fixed; \ 

w herein said method doesViot involve assembly of two or more fragments to form the ligand. 
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(Amended) 14. \ The system of claim 1 1 , wherein said means for performing the pre- 
docking conformational search comprises: 

means for randomly generating a plurality of conformations of the ligand; 
means fo\minimizing a strain of each conformation of the plurality of 
conformations; 

means for dsi ng the strain and a solvent accessible surface area of each 
conformation to rank the conformations; and 

means for clustering the conformations and retaining a desired number of top 
clusters of conformations, said retained number of top clusters of conformations 
comprising said multiple solution conformations of the ligand. 



(Amended) 16. The system of claim 15, wherein said means for generating the binding 
site image further comprises: 

means for placing a grid around the binding site of the protein; 

means for determining a hot spot search volume using said grid; 

means for determining hot spots using a grid-like search of the hot spot search 
volume; and 

for each type of hot spot, means for clustering the hot spots and for retaining a 
desired number of top clusters of hot spots, said desired number of top clusters 
comprising said multiple hot spots to be employed by said matching. 



(Amended) 18. Yhe system of claim 17, wherein said determining the unique rigid body 
K transformation comprises determining the unique rigid body transformation that minimizes: 




I(R, T) = £ L - RAj 
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I(r\j) = rms deviation between a j th hot spot and a j th atom of the at least one 
solution conformation; 

Hi \ position of a j th hot spot of the protein; 

position of a j th atom of the at least one solution conformation; 
a 3x3 rotation matrix; and 
a translation vector. 

(Amended) 1 9. I The system of claim 1 1 , wherein said means for optimizing comprises 
means for optimizing multiple protein-ligand complex formations, said means for optimizing 
comprising: 

mearfs for eliminating each ligand position having a predetermined percentage 
of ligand atoms with a steric clash; 

meanafor ranking remaining ligand positions using an atom pairwise score with 
a desired atom score cutoff, said atom pairwise score comprising a hydrogen bonding 
potential score \or a steric potential score; 

after ranking, means for clustering the ligand positions and selecting a top 
number n of ligand positions; and 

means fonoptimizing each ligand position of the n positions, allowing the 
translation, rotation and rotatable bonds of the ligand to vary. 



(Amended) 20. The system of claim 19, wherein said means for optimizing comprises 
means for optimizing each liaand position of the n positions using a Broyden-Fletcher- 
Goldfarb-Shanno (BFGS) optimization algorithm with a simple atom pairwise score, allowing 
the translation, rotation and ro\atable bonds of the ligand to vary. 



(Amended) 21 . At least one* program storage device readable by a machine, tangibly 
embodying at least one program obstructions executable by the machine to perform a 
method of docking a ligand to a protein, comprising: 

performing a pre-dockinp conformational search to generate multiple solution 
conformations of the ligand; 

generating a binding site \\age of the protein, said binding site image 
comprising multiple hot spots; 

matching hot spots of the binding site image to atoms in at least one solution 
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conformation of the Viultiple solution conformations of the ligand to obtain at least one 
ligand position relative to the protein; and 

optimizing the at least one ligand position while allowing translation, orientation 
A and rotatable bonds of the ligand to vary, and while holding the protein fixed ; 

wherein said method does not involve assembly of two or more fragments to form the ligand. 



(Amended) 24. nhe at least one program storage device of claim 21 , wherein said 

DCk\ 



performing the pre-dockfng conformational search comprises: 

randomly generating a plurality of conformations of the ligand; 
minimizing aWain and a solvent accessible surface area of each conformation 
of the plurality of confirmations; 

using the strainVrf each conformation to rank the conformations; and 
clustering the conformations and retaining a desired number of top clusters of 
conformations, said retailed number of top clusters of conformations comprising said 
multiple solution conformations of the ligand. 



[\ (Amended) 26. The at least one program storage device of claim 25, wherein said 

generating the binding site image further comprises: 

placing a grid around the binding site of the protein; 
determining a hot spot search volume using said grid; 
determining hot spots using a grid-like search of the hot spot search volume; 

and 

for each type of hot spot, clustering the hot spots and retaining a desired 
number of top clusters of hot spots, said desired number of top clusters comprising said 
multiple hot spots to be employed by said matching. 




^ / (Amended) 28. \ The at least one program storage device of claim 27, wherein said 
A k\determining the uniqu\ rigid body transformation comprises determining the unique rigid body 
transformation that minimizes: 
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I(R, T) = £ L - RAj - T 
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where: 

l(Rj) =y m s deviation between a j th hot spot and a j th atom of the at least one 

solution conformation; 

Hj = position of a j th hot spot of the protein; 
Aj = position of a j th atom of the at least one solution conformation; 
R = a 3*i3 rotation matrix; and 
T = a translation vector. 

(Amended) 29. The at least one program storage device of claim 21 , wherein said 
optimizing comprises optimizing multiple protein-ligand complex formations, said optimizing 
comprising: 

eliminating each\igand position having a predetermined percentage of ligand 
atoms with a steric c\ash\ 

ranking remaining^ligand positions using an atom pairwise score with a desired 
atom score cutoff,said atom pairwise score comprising a hydrogen bonding potential 
score or a steric potential score; 

after ranking, clustering the ligand positions and selecting a top number n of 
ligand positions; and 

optimizing each ligand position of the n positions, allowing the translation, 
rotation and rotatable bonds oftthe ligand to vary. 

(Amended) 30. The at least one program storage device of claim 29, wherein said 
optimizing comprises optimizing each ligand position of the n positions using a Broyden- 
Fletcher-Goldfarb-Shanno (BFGS) optimization algorithm with a simple atom pairwise score, 
allowing the translation, rotation and rotatable bonds of the ligand to vary. 

Attached hereto is a marked-up version of the changes made to claims 1, 4, 6, 8-10, 
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